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\ Vortex configurations in the two-dimensional torus are considered in noncommu- 

O , 

VO . solutions are constructed for the self-dual and anti-self dual cases by extending an 



o 



X 



Abstract 



algorithm originally developed for ordinary commutative space. We work within the 
Fock space approach to noncommutative theories and the Moyal-Weyl connection is 



Qh! used in the final stage to express the solutions in configuration space. 



1 Introduction 

The study of non-trivial classical solutions of field theories defined in noncommutative (NC) 
space-time has attracted much attention during the last years. The case of instantons, 
vortices and monopoles has been analyzed in great detail pQ, [2j. The analysis of these 
type of configurations simplifies for particular relations of coupling constants for which it is 
possible to establish the existence of Bogomolnyi-Prassad-Sommerfeld (BPS) equations |3j. 

For the particular case of vortices in Abelian Higgs model in the noncommutative plane, 
the existence of BPS equations was first established in As in the commutative space 

counterpart [Zj, no explicit analytical solutions exist. While in ordinary space the field 
profiles are found by solving numerically non linear differential equations [Zj, in NC space 
this is done by solving numerically non linear recurrence relations 0. 
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More recently, BPS equations for the Abelian Higgs Model in a two dimensional torus in 
NC space were found in jH]. As in the case of commuting coordinates, the analysis of solutions 
to these equations is more complicated than for the planar case due to the presence of non 
trivial boundary conditions. A very efficient numerical method was recently introduced in 
ordinary space by Gonzalez- Arroyo and Ramos jTT] . Guided by this method, we will address 
in this paper the problem of constructing numerical solutions of the BPS equations in the 
NC torus. 

This paper is organized as follows. In section 2 we introduce gauge and scalar fields 
defined in the noncommutative two-torus, discuss their boundary conditions, and analyze 
how gauge covariant and gauge invariant objects constructed from these fields should be 
integrated. The noncommutative Maxwell-Higgs model is introduced in section 3, where the 
derivation BPS equations and a Bogomolnyi bound for the energy are recalled Section 
4 is devoted to the construction of explicit vortex solutions to the BPS equations, this 
being achieved by using both the Moyal and the Fock space approaches for the treatment of 
noncommutative systems. 

2 Fields in the noncommutative torus 

We consider noncommutative 2 + 1 dimensional space-time with coordinates satisfying 

[x,y]=i6 (1) 



[x,t] = [y,t} = , (2) 

and the space coordinates are defined on a torus, (x, y) C T, the periods of T being (Li, L 2 ). 
We shall be interested in a U(l) gauge theory with a Higgs scalar <3> coupled to gauge fields 
Aj. The fields transform under the U(l) gauge group according to 

Ai -> AY = V~ l A i V + -V- 1 d i V i = l,2 (3) 

9 

<£ _> $ v = V' 1 $ , (4) 

with V G U(l) and g the gauge coupling constant. The fields are functions of (x, y), that 
is, they are operators acting on the Fock space generated by Eq.([T]). As we will be looking 
for static configurations, the time t, which in our approach is just a parameter, will not play 
any role. For definiteness, we will consider scalar fields in the fundamental representation. 
The other cases (anti-fundamental, adjoint) can be dealt in a similar way. Here, derivatives 
are defined as in the noncommutative plane, 

9i = } • (5) 
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As for the ordinary torus, a scalar field on the noncommutative torus will be defined as 
a function $(x,y) which is periodic up to gauge transformations. That is, 

$(x + L u y) = Ui{x,y)${x,y) = ¥ u ^\x,y) 

§(x 1 y + L 2 ) = U 2 (x,y)§(x,y) = ¥ u * 1 \x,y) , (6) 

where U\,U 2 C U(l) are the transition functions. Accordingly, boundary conditions for 
gauge fields are 

A l (x + L l ,y)=Af"\x,y) (7) 

Mx,y + L 2 )=Af* 1 \x,y) . (8) 

Consistency of the precedent relations leads to an equation for the U's which is, formally, 
the same as for the commutative torus, 

U 2 (x + L 1 ,y) U^x, y) = U x {x, y + L 2 ) U 2 (x, y) . (9) 

A solution of this consistency equations is given by 

E>i(x, y) = e l ™ L ^, U 2 {x, y) = e^^^ £ , (10) 

where 

u = —CL-s), keZ 
On 



s = V 1 - 2it6k/L 1 L 2 . (11) 

In the 9 — > limit, Eqs.(|10j l -()li p go smoothly to the solution corresponding to the commu- 
tative torus. 

Calling Ag the space of functions defined on T, a generic periodic function f(x, y) E Ag 
can be written in the form 



f&y) = fnm 2 (ni,n 2 ) , (12) 



ni,n 2 

where we have introduced 



(ni, n 2 ) — exp \2irin\— j exp ^iri^-jr-^J . (13) 

An integration in Ag, which we shall denote as Tr, can be formally introduced, 

I[f] = Trrf(x,y) = f 00 L 1 L 2 . (14) 
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Gauge invariant local objects are periodic in T and integrals of this kind of objects can be 
calculated according to this rule. Nevertheless, gauge covariant quantities f c , satisfy 

f c {x + L x , y) = U^x, y) / c (x, y) U^(x, y) 

f c {x,y + La) = U 2 (x,y) f c {x,y) U^(x,y) , (15) 

and it is simple to show that quantities such as f c are periodic in the scaled torus 7" with 
periods 

Li = sLi . (16) 

In this case the functions should be expanded in the basis 

x \ / y 
2mn->— 



(ni,n 2 )* = exp y2niriij- j exp y2iriri2j- J , (17) 

and the integral should be understood as 

I{f c ]=Tr f nx,y) = fc L 1 L 2 . (18) 

In theories defined in NC space, the more natural "local" (i.e, before integration) variables 
are covariant quantities (for instance the electromagnetic tensor, the energy density, etc). 
We will see then that the scaled torus T plays a fundamental role. 

Notice that the trace operation satisfies Tr (fg) = Tr(gf) and reduces to the standard 
integral on T in the commutative limit. One can see that this definition is crucial for 
preserving the cyclic property of the integral (trace) which in turn is essential in order to 
derive the equations of motion. For example, given two functions $i(x) and $2(2^) satisfying 
boundary conditions 

Qiix + L^y) = Uifay) ®i(x,y) 

$i(x,y + L 2 ) = U 2 (x,y)$i(x,y) , * = 1, 2 , (19) 

one can see that the product 

•11(f) $ 2 (f) (20) 
is strictly periodic in the torus T, but the transposed product 

$ 2 (f)$l(f), (21) 

is not periodic in T but in the scaled torus T. Nonetheless, as proved in [Hj, the cyclic 
property of the integral is still valid provided one integrates the first function in T while the 
second one in T 

Tit £2$) = Trr (^O?) H&) ■ (22) 
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That is, the cyclic property is preserved with the above definition. In what follows we shall 
indistinctly denote the trace operation by Tr assuming that the integrand is expanded in its 
natural domain of periodicity. 

So far, we have identified the space coordinate algebra defined in Eq. (0) with the algebra 
of creation-annihilation operators in a Fock space, and we have taken fields as operators 
$ in such a Fock space. As in the noncommutative plane, instead of working with fields 
depending on noncommuting coordinates Xi, one can work with ordinary coordinates Xi and 
introduce a noncommutative * Moyal product 

* x(x) = $(x) exp f^ u dX\ X (x) . (23) 

The connection between these two formalisms is found via the Weyl connection, an 
isomorphism that relates the algebra of functions multiplied with the noncommutative Moyal 
product and the algebra of operators in Fock space. For (x, y) G R 2 the relation reads 

$(x,y)4f(x,y) = $7v(x,y) , (24) 

where $(/ci,/c2) is the Fourier transformed of field <$>(x,y) defined in ordinary space. This 
formula can be easily extended to the torus. Indeed, as we signaled above, any function $ 
which is periodic in a torus T can be Fourier expanded as 

$(x,y) = ^2 ^nm 2 exp ^ 27rin i^J ex P \ 2nin 2 jr-J ■ (25) 



Then, Eq. (J24j) valid for R 2 , is replaced in the torus T by 

$(x,y) = }^ $mn 2 exp f — — 1 (m,n 2 > . (26) 

The connection between integration in both approaches is 

J d 2 x $(x,y) -> 2n9 Tr($) . (27) 

The Moyal mapping gives us the possibility to work with commuting coordinates. The 
difficulty with this approach is that the resulting expressions (and the equations of motion) 
are highly non local quantities in the sense that they involve derivatives of arbitrary order. 
As in the case of R 2 we will find more convenient to solve the equations of motion in the 
Fock space formalism, and use the Moyal correspondence to represent graphically the final 
results, by connecting operators with functions defined on configuration space. 
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3 The Maxwell- Higgs model 



We are interested in a model with a U(l) gauge field coupled to a Higgs scalar denned on 
the noncommutative torus. Dynamics of the model is governed by the Lagrangian density 

L=-\ + (A^) t " A ($ f $ - $o) 2 • (28) 

We are looking for static solutions to the equations of motion and hence we can look for 
minima of the energy (per unit length) 



E = Tr Q 4-4- + (A.$) t (A$) + A ($t$ - $2)2^) . 



(29) 



Here, 

= 4$ - igA $ (30) 
is the covariant derivative, and F%j is the electromagnetic tensor 

Fij = diAj - d j A l - ig[Ai, Aj] . (31) 

Notice that the covariant derivative used in Eq. (j3T7|) corresponds to a Higgs-gauge coupling 
which can be considered as in the fundamental representation (other choices are possible). 

As in the commutative case, the energy can be rewritten using the Bogomolnyi trick 
as |H], 

E =Tr Q |A$ - ijeij D^\ 2 + \ (4 -19£ij{® & - ^o))' + 

^A - y^j ($ f $ - $ ) 2 - 7 | $0 £ ij F ij + total derivative 

(32) 

where 7 = ±1, and E\2 = 1. 

Then, choosing coupling constants so that A = g 2 /2 (the Bogomolnyi point) one can 
establish a BPS bound for the energy 

£ > -7§<&oTr^4 , (33) 

where we have indicated that the trace is taken in the scaled torus f to stress that F^ is a 
covariant object. The bound is attained when the following BPS first order equations hold 

= 1 ge ij {®&- $1) . (34) 
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The sign of 7 should be chosen in such a way that the bound is positive. For definiteness, 
we shall set from here on 7 = — 1. 

In order solve these equations, let us start by observing that boundary conditions (jHJ) 
together with our choice of transitions functions (jlUj) imply for the gauge field the following 
relations 

Ai(x + Li,y) = Ai(x, y) 

Ax(x, y + L 2 ) = Ai(x, y) - ^thjL 2 

A 2 (x + Li, y) = A 2 (x, y) + -ixujLi 

9 

A 2 {x,y + L 2 ) = A 2 {x,y) . (35) 

A solution can be written as 

Ai(x,y) = Ai(x,y) + a { (x,y) , (36) 
where Aj are some periodic functions in the scaled torus T, and a, chosen as 

di = f Eiji 3 (37) 

with 

f = 

90 V *, 



f = -J^-~) ■ (3e 



The field strength Fj, can be written as 



^ij Fij ~\~ fij j (39) 



s 



where 



and 



27ik 1 

fij = en — yt~ ( 4 °) 
y Li\L 2 



F %j = diAj - djAi - i g [Ai, Aj] . (41) 
Here, we have introduced a scaled charge 

~g = sg. (42) 
Let us now study parameterization (|3fij) in connection with gauge transformations, 

AY = v- 1 [Ai + ai) V + - V~ l dt V 



V~ l AiV+ f £ij V' 1 x j V + - V- 1 di V . (43) 

9 



7 



Using Eq.flnj we can rewrite the middle term as a derivative term plus a, 

Aj = V- 1 Ai V - i Of V- 1 diV + di + - v- 1 di V 

9 



(44) 



= V- 1 Ai V + - V~ l diV + ai . (45) 
9 

Thus a gauge transformation on A\ is equivalent to a gauge transformation on A\ (keeping 
hi untransformed) but using the scaled charge g . 

We can summarize these results by stating that, in the analysis of gauge theories in the 
torus, one can trade non-trivial boundary conditions in the noncommutative torus T by 
periodic boundary conditions and a scaled charge g in a scaled noncommutative torus T. 

Let us now discuss the boundary condition equations for scalar fields. A field y) 
satisfying the boundary conditions © with transition functions given by Eq.QlOp. can be 
decomposed as 

$(x,y) = M~ 1 (x,y)x(x,y) , (46) 

where x(x,y) will be an explicit function fixed so that it satisfies the same boundary con- 
ditions as and M~ 1 (x,y), which has periodic boundary conditions on the torus T, 
will be found by solving the equations of motion. 
Thus x, must satisfy the conditions 

x(x +Li,y) = Ui(x, y) %(x, y) 

X(x,y + L 2 ) = U 2 (x,y)x(x,y) ■ (47) 

Using complex variables 

z = x + iy, z = x — iy , (48) 

a solution can be written as jH] 

\k 

X %1) =Ne%{*'*-*} Xldz^iz + aJ/I^iLx/L*) . (49) 

n=l 

Here 9^{z\t) is the Riemann 6 function 

O z (z\t) = J2e inrn2+2inB , (50) 



the ai are \k\ complex constants (which will be associated with the center of the vortices) 
satisfying 

1*1 

X>«=°> ( 51 ) 



n=l 
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and 

a = -— log(l -iru9) = -— logs . (52) 

In the 9 — > limit, this function coincides with the one obtained in the commutative case 
(see [TTj). Note that in the particular case of k = 1, Eq. ()49j1 simplifies to 

x(5,I)=AT e f{ i ' 5 - S }^(^/L 1 KL 1 /L 2 ) . (53) 



Using decomposition (|36| the BPS equations can be rewritten in the form 



2tt£; 



$$t 



7tu; 



where 



A z , Az 



(54) 
(55) 

(56) 



and the complex gauge fields are defined as A z = 1 2 2 ■ Since the fields are periodic in 
the scaled torus T, the total flux T of Fij on T vanishes (see equation (|14j) ). and then we 
have 



JF = Tr^ F 12 = Tr f f u 



2nk 



(57) 



It is easy to see that the ansatz ()46|) automatically satisfies the BPS equation ()55|) pro- 
vided that the gauge field is chosen as 



Az = -M l d s M 
9 



(5* 



where M is a (non unitary) function periodic in T. 

Then, it only remains to find M and M appearing in Eqs.(|49 p -()58j l so that the gauge and 
scalar fields satisfy the BPS equation (|54|). Defining 



// = MM f , 

the field strength F zz - can be written as 

F& = % -M- x H d z (H-%H^ M t_1 

and Eq . (|54j) takes the form 



Hd z (H- l d z H) = l -f [eH 



XX 



(59) 



(60) 



(61) 



where 

Ink 



6 = H - -=f±- ■ (62) 



g 2 L 1 L 2 

It is convenient at this point to introduce dimensionless fields and variables defined as 

$ — A, -> -^i— A u x % -> y/2g^i , (63) 



and redefine the parameters 
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A — > —A, 0^2<?& o e. (64) 

With these conventions, the Bogomonlyi point corresponds to A = 1, and Eq.(JHT|) becomes 

2 

H 4 (H- l d~ z H) = S -(eH-x X ] ) , (65) 
where we have redefined £ as 

47rA; , . 

e = 1 - . (66) 
L1L2 

In what follows we shall discuss in detail the numerical method used to solve this equation. 



4 Constructing solutions 

In order to find solutions to the BPS equations we shall first find H satisfying Eq. (j65|) and 
determine TV, the Higgs field normalization constant. We first need to compute, within the 
operator approach, the Fourier expansion of xic — ■N" 2 vtf > 

fjrjl = ^ Vn 1 ,n 2 (ni,n 2 )* . (67) 

ni,ri2 

Notice that the domain of periodicity is T, then the appropriate basis is 

(ni,n 2 )* = exp \ 2'nin x — ) exp f 2mn 2 — ) ■ (68) 
V sLxJ \ sL 2 J 

Surprisingly, it is possible to find as in the commutative space case, a closed expression for 
this quantity. Indeed, using the definition of 77 (Eq.([53|) for the single- vortex case k = 1), 
and after a long calculation, the coefficients ^„ 1>n2 (properly normalized as Trf (fjr)') = L\L 2 ) 
can be written as 



%i,n 2 = (-l) nin2 e 2s 'V ^ ; e 2 " 2min2 ^^. (69) 

To construct solutions to Eq.(|63)l we shall extend the technique described in jTTj for the 
commutative torus to the noncommutative case. Since for e = there is a trivial solution 

H = Constant, M = , (70) 
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we can use e as a perturbative parameter and expand H and the normalization constant M 
in powers of e 

oo oo oo 

H = Y,H k e k , H-^Y^H^, N 2 = Y,A k e k . (71) 

k=0 k=0 k=0 

Coefficients H k and H k , are operators periodic in T and can then be Fourier expanded, 

Hk= E^£n>^>*' = Y.~ h nlnM^T ■ ( 72 ) 

ni,ri2 ni,ri2 

Inserting these expansions in Eq. (j65|) one can determine order by order the coefficients, 



ni,na '"ni,n 2 



1 n\ — ri2 — 
ni^0,n 2 ^0 



(1) i m= n 2 = 

m,n 2 - i 2 ^, W2 i ^ ,n 2 ^0 

»?0,0 Knx.n^r ' ' 



where 



£m,n 2 = ttVT ( -^j= + mi ) , (74) 

with T = L2/L1 the aspect ratio of the scaled torus. In the same way one can calculate 
coefficients to any order N in e 

f rii = 122 = 

hiVno — S r {A) -r (B) -r {C) > (75) 

1, 2 I °"i."2 °"i."2 ni -^o,Ti 2 ^0 

with 

AT-l 

E 

ni,n 2 fc=l 

AT-l 

E 

m,n 2 fc=0 



^<3l,g 2 ~~ ^ ] ^ ] ^ l ni,n 2 ^'qi-n 1 ,q2~n2^Qi,Q2^qi~ni,q2-n2 ex P ( ^ ~ ~ ^((Zl ^l) 
ni,n 2 k=l x 1 z 

-kf]qi-n 1 ,q 2 -n 2 ex P I ^^2(^1 — n l) ) 

ni,n 2 k=0 " V ^1^2 / 

AT-2 y 47T 2 i9 

^91,92 = ^ ] ^ ] ^ l ni,n2^ l qi-n 1 ,q 2 -n 2 ^qi,q2^qi-n 1 ,q2-n2 ex P f ~~ ^((Zl — ^l) 

7 ^ \ -^1-^2 

ni,n 2 fe=0 1 z 

Coefficients /i nii „ 2 , appearing in the expansion of are obtained from HH^ 1 = 1, 

AT 

CI = - E E h nlnJ^'nL-n, ^ (^Mqi - n 1 )9 / L^) . (76) 
ni,n 2 fc=l 
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One also has to find a recurrence relation for the coefficients appearing in the expan- 
sion of the Higgs field normalization (Eq.(|71j)). For this, one uses the condition Tr f (-F12) = 
finding 

A = 

1 



A 1 



Vo,o 



An = V V V-m -mKi \n2 A N-k exp -- jr-~-ni7i2 , AT > 1 . 

(77) 



ni,n2 fc=l 



All these equations can be solved recursively. Assuming M Hermitian, it can be expanded 
in powers of e as well as in Fourier modes, and it is possible from H = M 2 to find recurrence 
relations for its coefficients. Finally, once the Fourier coefficients of H, H -1 , M, M~ l are 
known, together with the normalization constant Af, we can use decompositions ()46|) and 
()58|) to calculate all the fields in Fock space. The Weyl connection for periodic functions 
(J2HJ) can then be used to establish the correspondence between operators 0(x,y) and their 
associated functions 0(x,y). 

Using Eq. (j54j) for 7 = — 1 and the fact that Tr^F 12 = 0, it follows that 

Trf ($1 - 2 ? k ~ ) = Tr?$& > . (78) 



g 2 L\L 

Then 



slhh > 2 -4 , (79) 

g 2 

or in terms of the dimensionless variables, 

A > 4nk ^1 + Q = A c . (80) 

Then, the area A of the torus T has to be larger than the critical value A c in order for 
solutions to exist. 

We will focus first on the single- vortex case k — 1, and then make a few comments on 
k > 1. For simplicity we will only consider squared torus (Lj = L 2 ). We show in Fig.([TJ 
the solution for 6 = 2 and A = 100. Being the solutions periodic in T, we have represented 
them as a lattice of 9 plaquettes being the vortex solution centered in each plaquette. We 
plot both quantities F12 and $ * $t (the functions associated through the Weyl connection 
to the operators F 12 and in Fock space). For comparison, we also show the solutions 
for the same area and 6 = 0. In both cases, F12 has a maximum at the center of the torus 
(the location of each vortex). Unlike the commutative case, is different form zero at 
that point. 
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Figure 1: We represent F\2 and $ * & for 6 = and 6 = 2, for a torus of area A = 100. We 
consider 3x3 unitary cells, this leading to an array of 9 vortices. The distance among vortices 
equals L, which explicitly depends on 9. The magnetic flux is always 2irk, so as in the commutative 
plane, when incrementing 6, the vortices change conserving this quantity. 
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One can study the dependence of the solutions with the area for a fixed 9. We show in 
Fig.®, F 12 for several areas A > A c = Ank(l + §) for a fixed 9 = 2. Notice that for A = 300 
the configuration is already similar to the result in the noncommutative plane (see Ref.j^j). 
Alternatively, we can fix the area A and study the behavior of the solutions for different 
9 < 9 C = ~ 4 f fc . We show in Fig.(|2j) the results for A = 100 and several values of 9 in the 
range < 9 < 9 C ~ 14. 

It is also possible to construct gauge invariant quantities related to the scalar and gauge 
fields. Consider for example 

$t<£ = ^tfH-^f) , (81) 
and the invariant magnetic field B (to be distinguished form the covariant one F 12 ) 

B = (!>+<£) _1 ($^12$) = (&F 12 $)(&$)- 1 = -(1 - . (82) 

In the 9 — > limit, both quantities reduce to their analogues in commutative space. The 
flux of B across the torus now depends on 9 

Tr T {B) = 2n(l + ~^k. (83) 

We show in Fig.®, B configurations in the same conditions of Fig. (JTJ) . As it is gauge 
invariant, it is defined on the same torus T for all values of 9 (see Fig.pjl). 

We can also consider negative values for the noncommutative parameter 9. As the equa- 
tions remain unchanged when the following quantities are redefined as 

7^-7 , x 2 — ► — x 2 , A 2 — > -A 2 ,9^-9 ,k —k , 

this is equivalent to study solutions of the anti-self dual equations (this is, BPS Eqs.(|M jl -(|33 |l 
with 7 = 1) but with positive 9 parameter. 
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Figure 3: We show B configurations defined on a torus with area A = 100 for 8 = (left) and 
8 = 2 (right). The main difference with F\% is that it is defined on T, and its magnetic flux depends 
on 8. 



O 


o 





Figure 4: Upper views of B on a torus with area A = 100, for different values of 9. When 
9 -> 9 C ~ 14, B -c i 
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Figure 5: In the left we represent a two- vortex configuration of F\i for a torus with area A = 300, 
both for = and 6 = 2. The distance between vortices was fixed to 0.3L. In the right we show 
upper views of B in the same conditions. 

In the noncommutative plane, it has been shown that there exists a critical value 6* such 
that for 9 < 6>* no solutions to the self-dual equations exist [5], [TU]. In the units used in this 
paper, in the planer case this corresponds to 9* = — 2. The question that arises is if such 9* 
exists also in the NC torus and if it depends on the area. 

We have analyzed this problem numerically and we could not make the method to con- 
verge for 9 < —2, irrespectively of the value of the area. This is completely analogous to 
what happens in the noncommutative plane [3], jH], [HI], indicating that 9* = —2 also for 
the torus. Thus, we have been able to find solutions of the self-dual equations 9* < 9 < 9 C . 
Incidentally, notice that for 9 = 9*, the critical area is zero. The case of anti-self dual 
equations can be considered using the transformation mentioned above. 

It is also possible to obtain solutions for k = 2 and higher. In such cases we do not 
have closed expressions for the Fourier coefficients of 7717^ , but numerical calculations are 
straightforward. As an example, we show in Fig.© -F12 and upper views of B for both 9 = 
and 9 = 2. 

Finally, let us say a few words about the efficiency of the method. In order to analyze 



16 



the accuracy of the algorithm, we have rewritten Eq.(|65|). as 



e s 



- H- x N*Trf + %d z [H-%H] =1 



(84) 



and verified that the fourier coefficients A, 



of the LHS of this expression 



.4, 



91,92 — 




M 2 h 




(85) 



satisfy 



< 10~ p 



(86) 



for a given p. This is attained by increasing the number of Fourier coefficients and the order 
of the perturbative expansion in e. The same was done for the HH^ 1 = 1 constrain, and for 
other relevant equations. 

Convergence is slower for large area A ~ 20A C because, as solutions are more locali- 
zed, more Fourier coefficients are needed, and besides, e ~ 1. For instance, in this case 
a convergence for p = 3 is achieved with 441 Fourier modes, and iV = 300 orders of the 
expansion (this demands about 72 hours in a standard PC). As the area is reduced, the time 
of computation is considerably lower, being of about 20 seconds for e ~ 0. 

As an independent test, we have also calculated the magnetic flux and energy using the 
solutions and compared the results with the analytical values finding an agreement better 



than (~ 10" 4 ). 
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